Reproducing Z c (3900) through the ISPE mechanism 
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Being stimulated by the recent BESIII observation of a charged charmonium-like structure Z t .(3900), in this 
work we study the distributions of the 7/i/ot* and 7t + n~ invariant mass spectra of the y(4260) — > n + 7t~J/i// 
decays by the initial single pion emission (ISPE) mechanism, where the interference effects of ISPE mechanism 
with two other decay modes are also taken into account. The obtained dT(7(4260) — > n + 7i~J/if/)/dnij/^± and 
aT(Y(4260) — » n + n~J/if/)/dm n + n - marvelously agree with the BESIII data to reproduce the Z c (3900) structure. 
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Very recently the BESIII Collaboration announced the ob- 
servation of a charged charmonium-like structure Z f (3900), 
which is near the DD* threshold and comes from the analy- 
sis of the J/if/^ invariant mass spectra in e + e~ — > n + iC Jj\\) at 
-\fs - 4260 MeV (Jij. Before this interesting and important ex- 
perimental observation, there have been some theoretical pre- 
dictions of charged charmonium-like structures near DD* and 
D*D* thresholds yfl], which were mentioned in the BESIII 
paper on Z C (3900)W- 

The BESIII's result on Z c (3900) has also inspired exten- 
sive discussions on its properties Il5l-ll2ll and the prediction of 
a charged charmonium-like structure with hidden-charm and 
open-strange 11311 . What is more important is that the Belle 
Collaboration confirmed Z e (3900) in e + e~ — > n^n^Jlifi at 
yfs = 4260 MeV Q- Later, the former CLEO-c group ana- 
lyzed the data of their own and observed the charged Z c (3900) 
structure in their analysis of 586 pb~' data taken by their de- 
tector at if/ {4160) IU5I1 . which is well consistent with the BE- 
SIII's observation (Tj]. In addition, the neutral Z c (3900) was 
also released at a 3<x significance IU5I1 . 

Among many theoretical papers [2-4] on the predictions 
of charged charmonium-like structures cited in BESIII's pa- 
per |Q]], Ref. yfl has given the predictions of enhance- 
ment structures near DD* and D*D* thresholds in the cor- 
responding J/tf/n + , if/(2S)n + and h c {\P)n + invariant mass 
spectra, where 1^(4040), (K4160), ^(4415), and 7(4260) de- 
cays into 7T + n~J/if/, 7T + n~i(/(2S), and n + n~h c {\P) via the ini- 
tial single pion emission (ISPE) mechanism llql . which 
was first proposed to understand the Belle's observations 
of two charged bottomonium-like structures Zft(10610) and 
Z/,( 10650) appearing in the hidden-bottom dipion decays of 
T(5S) lfl7ll . Later, the ISPE mechanism was applied to 
study the hidden-bottom dipion decays of T(11020), where 
the charged bottomonium-like structures near the BB* and 
B*B* thresholds were predicted IU8ll . In addition, two charged 
strangeonium-like structures observable in the 7(2175) — > 



7r + 7r ^(1020) process were predicted in Ref. fll9ll using the 
ISPE. These abundant theoretical results predicted by the 
ISPE mechanism Jj, [jjl U9I1 can be further tested in experi- 
ment. 

Having more experimental information on the distribution 
of the 7/(/'7r ± and n + 7T~ invariant mass spectra of 7(4260) — > 
n + 7T~J/if/ (II, [3. U5I1 . we can perform an in-depth study of 
7(4260) — > n + 7T~J/if/ through the ISPE mechanism to get 
deeper understanding of a Z c (3900) structure discovered in the 
J/il/n* invariant mass spectrum. Along this road, in this paper 
we will fit our model based on the ISPE mechanism to the 
experimental data of the J/tfm* and 7r + 7r~ invariant mass spec- 
tra of 7(4260) — > 7t + 7t~J/i[/. In this study, we would like to 
check whether Z c (3900) structure can be reproduced through 
the ISPE mechanism. This will be helpful to deeply under- 
stand how to generate Z c (3900), which makes our study pre- 
sented here an intriguing research work. 

As for the hidden-charm dipion decay discussed here, 

7(4260)(p ) -» n + ( Pl )n-(p2)JmP3l 

there exist different decay mechanisms as shown in Fig. \T\ 
which play an important role to connect the initial state 
7(4260) and finial state 7r + 7T Jjifj. Among these decay mech- 
anisms, the first one is due to the 7(4260) direct decay into 
tt+tT J/if/ (see Fig. \T\ (a)), where 7(4260) -> 7t + n-J/tf/ oc- 
curs without the intermediate state contribution and mainly 
provides the background for the distributions of the ti + ji~ and 
J/ifrn* invariant mass spectra. The intermediate state contri- 
butions come from cr(600), /o(980), and the triangle hadronic 
loops composed of charmed mesons. 

By the effective Lagrangian approach, the decay amplitude 
of the direct decay of 7(4260) — » n + n~Jlif/ can be written as 



M Dire ct[ 7(4260) -> n^-J/ilj} 
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which was once constructed by Novikov and Shifman to study 
if/' — > 7t + n~J/tf/ transition 112011 . In Eq, ([T), the subscripts, 
S-wave and D-wave, represent the corresponding amplitudes 
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FIG. 1 : (Color online.) The typical diagrams depicting 7(4260) — > J/ipn + n~ decay. Here, Fig. Q] (a) denotes 7(4260) direct decay into 
,//i/OT + 7r~; Fig. [3(b) describes the intermediate hadronic loop contribution to 7(4260) — > J/i(/7r + n~. Fig. [3(c) and (d) are from the ISPE 
mechanism. 



from S-wave and D-wave contributions, respectively. AM — 
mY(4260) - tnj/ij, is the mass difference between 7(4260) and 
J l\f/. q 2 = (p\ + p 2 ) 2 = m\+ n - denotes the invariant mass of 
;r + 7r~ while 9 is the angle between 7(4260) and ii~ in the n + 7t~ 
rest frame. In addition, the pion mass and decay constants are 
taken as m„± = 139 MeV and f„ - 130 MeV, respectively. 
In the direct decay process, two free parameters, F and k, are 
introduced, which can be determined by fitting to the experi- 
mental data. 

The second mechanism in the 7(4260) — * Jlipn + n~ decay 
is due to intermediate S = <x(600) and /o(980) contributions, 
where the hadronic loop constructed by the charmed mesons 
is the bridge to connect 7(4260) with j/ipS as described by 
Fig. [T](b). In this work, we parameterize the decay amplitude 
corresponding to Fig. [T|(b) as 



M S [Y(426Q) -> tt + 7t-J/^] 



= 4 



F(4260) e y/>A Sfy 



Fs 



(Pi + Pi) 1 -m\ + im s V s ' 



(2) 



where a subscript >S can be taken as cr or /o(980) and ms and 
Y$ are the resonance parameters of cr(600) and /o(980). Since 
the total decay width and mass of <x(600) are of the same or- 
der, the minimal width approximation is not valid. Here, we 
adopt the momentum dependent decay width in the propaga- 
tor of cr(600), i.e., 



^a-{m n * n -) = T 



m a \p(m^ jr )\ 
m^jr \p{m a )\ 



(3) 



where \p(m n + n -)\ - Jm^ r /4 - m 2 is the value of pion mo- 
mentum and \p(m cr )\ is the pion momentum with on-shell cr 
meson. In Eq. ©, two fitting parameters F a and F '/ (980) are 
introduced. 

Besides these two decay mechanisms mentioned above, 
there is the third mechanism of 7(4260) — > //</ot + 7t~ decay, 
i.e., the ISPE mechanism, via which the typical diagrams de- 
picting 7(4260) — > J/tpn + n~ decay are given in Figs. Q](c) and 
(d) corresponding to the n~ and n + emissions from 7(4260), 
respectively. Through the ISPE mechanism, 7(4260) first de- 
cays into a pair of intermediate charmed and anti-charmed 
mesons with initial single pion emission. Due to the continu- 
ous energy distribution of this pion, the intermediate charmed 
and anti-charmed meson pair is of low momentum, which 
makes them easily interact with each other to transform them 
into J /if/ and n. 



The detailed calculations of 7(4260) — > n + n~J/ip via the 
ISPE mechanism are presented in Ref. yfl, which shows that 
there are charged enhancement structures near DD* and D*D* 
thresholds. After performing the loop integrals, the decay am- 
plitudes via the intermediate DD* + h.c. and D*D* can be pa- 
rameterized as 

Mf s f E [7(4260) -> JT + 7T-J/i//] 

- gY(426Q)D'DK^r {4 2(,Q^ V j/,i,[ A 0gfiv + (#lPl/i/?lv 

+B 2 p\nP2v + B 3 p 2lt p\ v + B 4 p2t l p2y)\, (4) 

Atgpf [7(4260) -> 7T + 7T-J/t/s] 

= gY(4260)D'D'7r^ ( 4 26O - ) ej /l// \A' g llv + [B'iPlfiPlv 

+B' 2 pi ll p 2v + B' 3 p 2ll pi v + B' 4 p 2lt p 2v )j, (5) 

respectively, where the coefficients Aq and B. (i = 1 ~ 4), in 
front of the different Lorentz structures can be evaluated by 
the loop integrals in the ISPE mechanism (see Ref. yfl for 
more details). As indicated in Ref. 131, the coupling con- 
stants of 7(4260) interacting with D ( *'D ( *V, gY(4 2 60)D'DK and 
gY(4 2 60)D'D'7i, are unknown and become fitting parameters. 

The total decay amplitude of 7(4260) — > n + n~J/ip is the 
sum over the subamplitues, i.e., 

M Total [7(4260) -> 7r + 7T-J/if/) 

= Aloirect + e^Mo- + e''^ 8 °>M /o{9 80) 



+e 



'0ISPE 



l^ISPE + ^ISPE j ' 



(6) 



where three phase angles (pa-, 0/ o (98O) and </>ispe are introduced. 
As discussed above, in our scenario we introduce nine free 
fitting parameters as 

F, K, Fa-, F/ (980), (pa; 0/ o (98O)> 0ISPE 
gY(4 2 60)D'Din gY(4 2 60)D*D*n, 

which are determined by fitting to the experimental data. 

With the above preparation, the differential decay width of 
7(4260) -> TrW/iA is 



dT (7(4260) 
1 1 



-> TT+TrV/l/O 

I . -TVrfol 1 2 



3 (2tt) 3 32ml 



M ToUl \ Z dm 2 



Jl*n dm nn> 



(7) 



7(4260) 



where m 2 J/i/m = (p 2 +Pi) 2 and m 2 ^ = {p\ +p 2 ) 2 . Using Eq. ©, 
we can obtain the theoretical distributions of the J/ipn* and 



7r + 7T~ invariant mass spectra, which will be applied to fitting 
to the BESIII data |Q]]. Additionally, the resonance parameters 
involving in our calculation are listed in Table H 



TABLE I: The resonance parameters (in units of GeV) used in 
present work. The resonance parameters of <r(600) arc taken from 
Ref. 12111 and the other resonance parameters are from PDG [22]. 

Parameter value Parameter value Parameter value 

»Jr(4260) 4.263 nij/tp 3.097 m„± 0.139 

m a 0.513 I"V 0.335 

m /o(9go) 0.980 r /0(9 g ) 0.070 

At present, BESIII has not only given the J/tfm* 1 invariant 
mass distributions but has also provided the n + 7T~ invariant 
mass spectrum |l||. This information provides us good chance 
to test the ISPE mechanism. With Eq. (|7]i and the help of 
the MINUIT package, we carry out the global fit to the BE- 
SIII data of 7(4260) -> n^n^J/ip, which is helpful to answer 
the question whether the Z c (3900) structure observed by BE- 
SIII can be reproduced by the ISPE mechanism, where we 
also consider the interference effects of the ISPE mechanism 
with other two decay mechanisms. We fit the <fF(F(4260) — > 
7T + 7r~J/i[f)/dmj/j,„± and dT(Y(426Q) — > n^'J/i^J/dm^^- dis- 
tributions with the experimental data. 



TABLE II: The obtained values of the fitting parameters for the best 
fit to the BESIII data of F(4260) -> TT+ir-J/if/ decay [J. 



Parameter 



Value 



Parameter 



Value 



24.30 ± 6.05 



K 



190.27 ± 28.35 GeV 2 



0.25 ±0.15 
-1.09 ±0.51 

/ A(980) 83.56 ±11.54 GeV 2 /o(98O) -2.52 ±0.23 

gYmmo-Dn 25.26 ±1.88 g m2 60)D'D^ -2.30 ± 1.11 GeV" 1 
1SPE 2.87 ± 0.30 

The values of the fitting parameters obtained are listed in 
Table [TT] The corresponding best fit to the the J/ifm* 1 and 
7r + 7r~ invariant mass distributions is shown in Fig. [2] Com- 
paring our theoretical results with the BESIII data, we no- 
tice that the oT( 7(4260) — > Ji + ji^JI^/)ldmji^„± distributions 
well agree with the J/if/n* invariant mass spectra. Here the 
Z c (3900) structure and its reflection can be nicely reproduced, 
where our fitting results indeed reflect the fact that Z c (3900) is 
much sharper than its reflection. In addition, with the same fit- 
ting parameters, the theoretical distribution of aT(7(4260) — > 
n*n~ J l\(/) I drrix*^ is also given, which can well reproduce the 
BESIII data of the tt + ji~ invariant mass spectra again. 

In summary, the charged charmonium-like structure 
Z c (3900) has been newly observed and reported by the BESIII 
Collaboration 01 and confirmed by Belle OH and CLEO-c 



IU5I1 . After this new observation, different theoretical groups 
have given different explanations to Z c (3900) |p|- tl2ll . Among 
these explanations, the discussions on whether the structure 
of Z c (3900) is due to an exotic state or not are very popu- 
lar, where the exotic state explanation of Z f (3900) includes 
molecular state composed of D and D* mesons or tetraquark 
state. 

As mentioned in the BESIII experimental paper [lj], there 
exist the theoretical predictions of charged charmonium-like 
structures near the DD* threshold 12-4-] before the observa- 
tion of Z f (3900). In Ref. Q], the ISPE mechanism was ap- 
plied to study 7(4260) — > n^n^J/iff, ji + n~\jj' , n + n~h c decays, 
where two of the authors of this paper predicted the interest- 
ing phenomena of charged DD* and D*D* structures exist- 
ing in the J/tj/n*, 4>' n± an d ^V* invariant mass spectra. The 
BESIII observation of Z f (3900) has inspired us to study the 
7(4260) — > n + n~JI4i decay process by the ISPE mechanism 
again together with other diagrams shown in Figs. Q] where 
we are able completely to fit our model including a tree dia- 
gram as well as relative phases to the experimental data in our 
scenario since BESIII now provides abundant information of 
the distribution of the J/ifm* and n + n~ invariant mass spectra. 

As shown in Fig. [2] the Z c (3900) structure has been remark- 
ably reproduced under the ISPE mechanism that includes the 
interference effects. This fact affirmatively answers whether 
Z c (3900) corresponds to the predicted charged charmonium- 
like structure near DD* threshold via the ISPE mechanism in 
Ref. H. 

Besides the prediction and/or reproduction relevant to 
Z c (3900) under the ISPE mechanism, there are other abundant 
predictions on the features of charged enhancement structures 
in the hidden-charm dipion decays of higher charmonia y], 
the hidden-bottom dipion decays of T( 11020) IU8I1 . and the 
hidden-strange dipion decays of 7(2175) IU9I1 . We hope that 
future experiments will carry out search for these novel phe- 
nomena. 

Given the experimental data, the present work has provided 
an important approach to test the ISPE mechanism. With more 
experimental observations, more phenomenological and theo- 
retical efforts should be made by combining further experi- 
mental data, which will finally clarify what is the underlying 
mechanism behind the observed Z c (3900) structure and other 
exotic structures. 
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FIG. 2: (color online). The distributions of the J/iftn + , J/ifm' and n + n~ invariant mass spectra of K(4260) — > n + n~JI>l/. The blue dots with error 
bars are the experimental data given by BESIII QO, while the red histograms are the best fit by our model considering the ISPE mechanism 
contribution to the K(4260) — > tc + 7x~JI>I/ decay. 
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